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As photoelectrically detected *’Cf-source-driven neutron signals always contain noise, a denoising algorithm 
is proposed based on compressive sensing for the noised neutron signal. In the algorithm, Empirical Mode 
Decomposition (EMD) is applied to decompose the noised neutron signal and then find out the noised Intrinsic 
Mode Function (IMF) automatically. Thus, we only need to use the basis pursuit denoising (BPDN) algorithm 
to denoise these IMFs. For this reason, the proposed algorithm can be called EMDCSDN (Empirical Mode 
Decomposition Compressive Sensing Denoising). In addition, five indicators are employed to evaluate the 
denoising effect. The results show that the EMDCSDN algorithm is more effective than the other denoising 
algorithms including BPDN. This study provides a new approach for signal denoising at the front-end. 
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I. INTRODUCTION 


The photoelectric detection and imaging technology have 
been widely used in Nuclear Material Identification System 
(NMIS) to detect neutron signal and nuclear components. 
Time-frequency method has been used to analyze neutron sig- 
nal and identify properties of nuclear material [1, 2], while to- 
mographic imaging has been used to provide more informa- 
tion about geometry of nuclear components, which leads to 
high accuracy of NMIS [3-5]. However, it is inevitable that 
there are noises in the measurement process due to external 
environment, detector and electronic device, and the noises 
may confuse some useful but weak information. Thus, we 
should denoise the neutron signals first to improve accuracy 
of NMIS. 

Except the circuit and technic neutron signal denoising 
methods, traditional denoising algorithms mostly represen- 
t the noised signal in transform domain and then threshold for 
the purpose of denoising. As a result, FFT filter and wavelet 
shrinkage has been applied in neutron signal denoising [6]. 
Besides, a new denoising algorithm named EMDSD (Empiri- 
cal Mode Decomposition Double Smoothness Detecting) has 
been proposed for neutron signal [7]. However, these denois- 
ing algorithms cannot adjust decomposition base automati- 
cally. 

In recent years, the paradigm of sparse sampling and recon- 
struction, called compressive sensing (CS), is a state-of-the- 
art research [8, 9]. CS approaches have opened up many new 
research avenues in the field of under-determined systems, 
and found many practical applications in image processing, 
wireless communication, data-streaming, and medical reso- 
nance. Also, CS has been applied in image denoising [10- 
14], and some new denoising algorithms [15, 16] are derived 
from basis pursuit denoising (BPDN) algorithm [17]. Al- 
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though these algorithms are proposed for all signal frequen- 
cy, the compressive reconstruction error will still increase and 
signal may be noised again. 


Fortunately, the Empirical Mode Decomposition (EMD) 
can decompose the signal into Intrinsic Mode Function (JM- 
F), from finer temporal scales (high frequency IMFs) to 
coarser ones (low frequency IMFs). As the noise is always 
in high frequency IMFs, BPDN algorithm can be adopted 
to denoise them. In this paper, we propose a new denoising 
algorithm named EMDCSDN (Empirical Mode Decomposi- 
tion Compressive Sensing Denoising). Several denoising al- 
gorithms are compared with EMDCSDN algorithm to verify 
its robustness. 


Il. METHODOLOGY 
A. Materials 


Neutron signals from *>*Cf-source-driven nuclear material 
fission are detected by photoelectric detectors. We designed 
a *>°Cf-source-driven verification system to identify proper- 
ties of 5U (Fig. 1(a)). It consists of a 252Cf neutron source, 
an ionization chamber, scintillation neutron detectors, work- 
station with high-speed (1 GHz) data acquisition card with a 
large-capacity disk array, verification software and user in- 
terface (UI) system. Three detectors are placed around the 
fissile material. The target-detector distance d and the angle 
a between detectors can be adjusted according to the mea- 
surement requirements (Fig. 1(b)). 


Figure 2 shows the noised neutron pulse signal acquired by 
the ***Cf-source-driven verification system. For experimen- 
t purpose, the EMDCSDN algorithm is adopted to denoise 
this signal. As the pure neutron signal is unknown, the curve 
fitting of noised neutron signal is regarded as the pure signal, 
so as to evaluate the denoising effect. All the materials are 
used without further purification. 
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Fig. 1. (Color online) Prototype of 252Cf-source-driven verification system (a) and position of detectors (b). 
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Fig. 2. (Color online) The noised neutron signal. 


B. Basis pursuit denoising 


Compressive sensing, the paradigm of sparse sampling and 
reconstruction, is an advanced theory. Specifically, let x= 
{x1,%2,...,¢n},x€ RY with x= Wa, where a has only K 
non-zero elements and K « N. The x is called as K-sparse 
with respect to the transform W. So the random measurements 
y= {y1, y2,---, yn}, ye R” are generated by 


y=@x, where 8e RYN, (1) 


where, ® is a randomly sampling matrix, and the number of 
measurement M < N. It is an ill-posed problem recovering 
x from y . In the compressive sensing theory, it testifies that 
K-sparse signal x can be recovered by M = O[K log(N/K)] 
measurements as long as ® satisfies the restricted isometry 
property (RIP), and the reconstruction can be achieved with 
probability close to one by solving the following convex op- 
timization. 


a=argmin|la||;, subject toy = Pya, (2) 


where ||a||ı denotes the /;-norm of the vector a. 

Thus, the compressive sensing process of neutron signal 
can be described as follows. It assumes that the noised neu- 
tron signal x= s +n = pa + n, where s is the original 


signal, n denotes additive Gaussian white noise and w is the 
sparse basis. Consequently, the compressive sampling can be 
defined as, 


y=@(s4+n)=6s+On=6s+2=GHya+z, (3) 


where z is the noise in sampling and ||z||2 < €. 

The reconstruction can be achieved by Eq. (2). With noisy 
or imperfect data, it is impossible to fit the linear system ex- 
actly. Instead, the constraint in BP is relaxed to obtain the 
basis pursuit denoising problem. 


min |la||; subject to ||y — ®wall2 < e. (4) 


Therefore, an effective algorithm using Spectral Projected 
Gradient (SPG) is used to solve BPDN [17]. It can solve the 
BPDN problem much more quickly, and the memory require- 
ments are constant throughout all the iterations. However, 
most of the noised signal frequencies are pure and there is 
no need to denoise them. So, if we use the BPDN algorithm 
for all frequencies, the signal will be noised again and some 
useful information may be filtered. Therefore, based on com- 
pressive sensing, a new denoising algorithm for the noised 
frequencies only is proposed. 


C. EMDCSDN algorithm 


Hilbert-Huang Transform (HHT) is a new method of 
time-frequency analysis technique, which was proposed by 
Huang [18], and the essence of this method is EMD. The 
EMD is advantageous in that the base functions are derived 
from the signal itself. Hence, the analysis is more adaptive, 
in contrast to the wavelet method where the base functions are 
fixed. Any signal can be decomposed into a finite and a small 
number of IMFs to give meaningful instantaneous frequency, 
as described by Eq. (5). 


N 


x= XO IMF; + res, (5) 
i=1 


where res means a residual of EMD. 
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Figure 3 shows the IMFs of neutron signal decomposed 
by EMD. It seems that only the first few high frequencies of 
IMFs are noised, and rest of the JMFs are pure. Although 
the EMD itself has characteristics of filter, the noised IMF's 
cannot be filtered directly. Therefore, we need to find the 
noised IMFs, to use BPDN algorithm for denoising. 
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Fig. 3. (Color online) IMFs of neutron signal decomposed by EMD. 


Thus, the EMDCSDN algorithm is proposed as shown in 
Fig. 4. 


Noised Denoised 


Fig. 4. The EMDCSDN algorithm. 


Algorithm 1: EMDCSDN algorithm for neutron signal 


Input: noised neutron signal x; 

Output: denoised neutron signal <p; 

Decompose x to IMFs IMF,, j = 1,2,..., N) and res; 
Select the noised IMFsUMF;,., c = 1, 2,...,C) and pure 
IMFsUIMF,, p = 1, 2,..., P) by IMF selector, where 
N=C+P; 

Use BPDN algorithm to denoise the noised IMFs, and the 
output denoised IMFs are DLIMF.,c = 1,2,...,C; 


Calculate rp = D D_IMF. + Dii IMF, + res; 


In this EMDCSDN algorithm, the IMF selector is designed 
according to the characteristics of white noise [19]. One of 
the characteristics is that I; = EP; is a constant, where F; 
is the energy density and P; is period of the IMFs. 


ie 2N 
B: = y LIC, P, = , (6) 


Nmax (i) + Mmin (i) 
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where, N, C, nmin and Nmax are the length, amplitude, the 
number of maxima and the number of minima of the IMFs, 
respectively. 

If R; = |T; —F)/P| > 1, then T}(j = 1,2,...,i — 1) 
is a constant and the first i — 1 IMFs are noised, where I is 
the mean of 7}. Thus, we can use this IMF selector to find 
out which JMF is noised. Figure 5 shows that R3 > 1, so 
for the noised neutron signal, IM Fı and IM F need to be 
denoised. Thereby, the EMDCSDN algorithm can narrow the 
denoising scope and reduce the compressive sensing recon- 
struction error, and the denoising effect will be improved. In 
the next section, how to evaluate the denoising effect will be 
discussed. 
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Fig. 5. (Color online) The JMF selector results of neutron signal. 


D. Denoising evaluation 


In general, the signal-to-noise ratio (SNR), peak signal-to- 
noise ratio (PSNR) and mean squared error (MSE) are adopt- 
ed to evaluate the denoising effect. However, only the three 
indicators are not enough because high SNR means that some 
useful information has been filtered. So, another two indica- 
tors, smoothness of curve (SOC) and correlation coefficient 
(CC) should also be involved. The SOC and CC can be 
defined as: 


N-1;~ A 2 
SOC = Liat BU EJ 1) ima a(t ] , (7) 
iar leli + 1) — (a)? 
N _ 
2 (wi — £)(4i — ĉ) 
cC= = (8) 
N N a 
Ee — 2) 2 (i — &)? 


where x is the original noised signal and ĉ is the denoising 
signal. For CC indicator, the bigger the better. However, for 
SOC indicator, smaller SOC means better denoising effect. 
Hence, in order to evaluate the denoising effect, we need to 
consider all the indicators above. 
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Il. RESULTS AND DISCUSSION 


For the noised neutron signal x, x € {R} with elements 
x[n], n = 1,2,...,N and N = 1024. The number of com- 
pressive sampling is regarded as M, then the corresponding 
compressive sampling rate is M/N. 

As mentioned before, the first two IMFs IMF}, IMF2) are 
noised. Therefore, BPDN algorithm is used to denoise the 
two IMFs. The results are shown in Fig. 6, where DN is the 
abbreviation of denoising and M/N = 0.65. 


o 100 200 300 400 500 600 700 800 900 1000 
T ¥ T T T T id T T * 
T 
=. 
m 
i 
= 
=z - 
a 
È 
í 
N 
= 
= 
z 
5 
1 n 1 n 1 1 1 L L n 1 n 1 n 1 
o 100 200 300 400 500 600 700 800 900 1000 
Time (ns) 


Fig. 6. (Color online) Denoised IMFs by EMDCSDN algorithm 
(M/N = 0.65). 


Figure 7 shows the EMDCSDN algorithm denoising result 
of neutron signal. However, this cannot fully reflect the de- 
noising effect, it should be evaluated by the five indicators of 
SNR, PSNR, MSE, SOC and CC. 

Different compressive sampling rates M/N shall lead to 
different denoising effects. Figure 8 shows the denoising 
effects of EMDCSDN algorithm for different M/N. Figure 8 
has double Y axis. The left represents the MSE, SOC and CC, 
and the right one is SNR. SNR and CC increase with M/N, 
while MSE and SOC decrease with increasing M/N. The 
best denoising effect can be seen at M/N = 0.75, where 
SNR and CC are the biggest, and MSE and SOC are the s- 
mallest, and the denoising effect is the best. However, if we 
want a lower compressive sampling rate, M/N = 0.65 may 
be a better choice because its denoising effect is very close to 
M/N = 0.75. Here, we choose M/N = 0.75. 

The EMDCSDN algorithm is compared with four denois- 
ing algorithms: (1) EMD Filter (EMDF) that filters the noised 
IMFs directly; (2) Basis Pursuit Denoising (BPDN), com- 
pressive sensing denoising problem solved by SPGLI algo- 
rithm; (3) Wavelet Threshold Denoising (WTDN) that de- 
noise signal in wavelet domain with threshold; (4) EMD 
Wavelet Denoising (EMDWTDN), which uses WTDN to de- 
noise the noised JMF’s and the reconstruct signal with all ZM- 
Fs. The denoising effects are given in Table 1. 

From Table 1, the EMDF algorithm is a simple but less 
effective algorithm because it also filters out some useful in- 
formation. As the EMDWTDN algorithm is better than the 
WTDN algorithm due to the EMD, we think that EMD can 
improve the denoising effect and the effect of EMDCSDN 
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Fig. 7. (Color online) Denoising result of neutron signal (M/N = 
0.65). 
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Fig. 8. (Color online) Denoising effect of EMDCSDN algorithm for 
different M/N. 


algorithm confirms this view. Considering all the five indi- 
cators, the EMDCSDN algorithm is more effective than the 
other algorithms due to the EMD and the noised IMFs auto- 
matic selector, although its processing time is a little longer. 
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TABLE 1. Contrast between EMDCSDN and another four denoising algorithms 


noised neutron signal 


Algorithm SNR PSNR 
24.24559 39.93027 
EMDF 18.16586 33.57995 
BPDN (M/N = 0.75) 27.59472 42.93629 
WTDN 28.74912 43.98175 
EMDWTDN 29.48255 44.98517 
EMDCSDN (M/N = 0.75) 30.62113 45.89118 


IV. CONCLUSION 


A new denoising algorithm, EMDCSDN, has been pro- 
posed for noised neutron signal. In fact, it is a modified 
BPDN algorithm because we only need to denoise the noised 
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